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Top quark studies at future e^e~ colliders are considered. Two issues are discussed: a - Some results are 
presented on the decays of top quarks. Energy distributions of charged leptons and neutrinos in t ^ bW he^v 
and jets in t ^ hW — > bdu decays are sensitive to the structure of thW vertex. Distributions of charged leptons 
from top decays are particularly useful in polarization studies whereas neutrinos are sensitive to deviations from 
the Standard Model, b - Recent calculations are reviewed on the top quark pair production in e^e~ annihilation. 
The differential cross sections in the threshold region can lead to an accurate determination of the top quark 
mass and the interquark potential. The effects of the top-Higgs Yukawa coupling and some higher order QCD 
corrections are also under control. 



1. TOP QUARK AT e+e COLLIDERS 

Top quark is the heaviest fermion of the Stan- 
dard Model. The existing direct lower mass 
limit||l| rrit > 131 GeV as well as the value 
rut = 174li[2ii9 GeV derived from the recent fit 
to LEP+SLC data |^ indicate that its production 
is out of reach of LEP II. (Two weeks after this 
talk had been delivered the top quark saga cul- 
minated at the press conference at FNAL when 
evidence for tt production was announced by the 
CDF Collaboration . The reported value of the 
top quark mass rrit — 174±10t^g GeV agrees very 
well with the indirect determination ||].) 

It has been argued that single top produc- 
tion may be observed for rrit up to about 180 
GeV if LEP II center-of-mass energy is pushed 
to -^/s = 210 GeV. However, a recent article [|j 
found a dramatic cancellation of contributions to 
the reaction e''"e~ — s- e~Dtb. A similar destruc- 
tive interference damps the single top production 
in 67 collisions [||. It is noteworthy, however, 
that this does not occur for pp and pp collisions 
|6| where for large top masses the cross section 
of single top production becomes comparable to 
that of tt production. In Fig.^ the cross sction 
of reaction e^e~ e~vth is shown for ^/s =170, 
190 and 210 GeV @. The solid fines correspond 
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to the complete set of Feynman diagrams, the 
dashed lines to the photon exchange subset and 
the dotted lines to the Weizsacker- Williams ap- 
proximation. It is evident that even for the in- 
tegrated luminosity C = 500 ph~^ and ^/s ~ 210 
GeV one expects no events. Therefore, the top 
quark can be observed in e^e~ collisions only if 
the Next Linear Collider is built. Comprehensive 
studies have been performed of the physics poten- 
tial of a 500 GeV center-of-mass energy collider. 
The results are summarized in . Theoretical 
studies on the top quark are reviewed in [^^-|l3| . 

2. TOP QUARK DECAY 

2.1. Total decay rate 

Top quark is the first heavy quark whose mass 
can be measured to better than 1% precision at 
a future e+e" collider |l^. Therefore, measure- 
ments of its width can test not only the Standard 
Model at the Born level, but also the QCD radia- 
tive corrections which are about -8.5%. This is in 
contrast to h and c quarks, where uncertainties in 
the masses and non-perturbative effects preclude 
this possibility. 

Being much heavier than the W boson, 
the top quark decays dominantly into bW fi- 
nal states. In the Standard Model with 
three quark-lepton families, which is assumed 
throughout, the modulus of the element Vtb of 



Table 1 



Top quark width as the function of nit including (Ft) and without (r„^) the contributions of the QCD, 
non-zero W boson width (WBW) and electroweak (EW) corrections. 



mt[GeV] 


rt[GeV] 




QCD[%] 


WBW[%] 


EW[%] 


150.0 


0.809 


0.885 


-8.47 


-1.69 


1.57 


160.0 


1.033 


1.130 


-8.49 


-1.60 


1.62 


170.0 


1.287 


1.405 


-8.49 


-1.52 


1.67 


180.0 


1.572 


1.714 


-8.48 


-1.45 


1.70 


190.0 


1.890 


2.059 


-8.47 


-1.39 


1.73 


200.0 


2.242 


2.440 


-8.46 


-1.33 


1.76 



Figure 1. Total cross sections for the reaction 
e+e^ — + e~ Dtb{e^ vib) as the function of the top 
quark mass |^. 



the Cabibbo-Kobayashi-Maskawa quark mixing 
matrix is close to 1 ( \Vtb\^ 0.9985 to 0.9995 ) 
. Neglecting small corrections due to the non- 
zero mass of b quark and the width of W boson 
the total decay rate can be written as follows: 
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The one loop electroweak corrections to the to- 
tal decay rate have been calculated in the nar- 
row width approximation p7[ | . They have turned 
out to be positive and rather small (1-2%). A 
reason for this is that the contribution from the 
Higgs-top quark Yukawa coupling remains very 
small for realistic top quark masses [|l^ . The ef- 
fect of the non-zero W width is comparable in 
size to the electroweak correction but of the op- 
posite sign jl^. Formulae for the QCD corrected 
total decay rate including non-zero b quark mass 
and W boson width are given in (l|Jl|l. In Table 
1^ we present the results jl^ for the total decay 
rate Ft and its narrow width Born approximation 
Fnu). We give also the Standard Model contribu- 
tions to the width of the top quark from the first 
order QCD corrections, W boson width (WBW) 
and the electroweak (EW) corrections. A number 
of intrinsic uncertainties remains. It should be 
noted that the size of the electroweak corrections 
is comparable to the uncertainties from the as yet 
uncalculated ©(a^^) correction. The present un- 
certainty in as and the ignorance concerning the 
second order QCD correction limit the accuracy 
of the prediction to about 1-2%. Experimental 
errors as well as theoretical uncertanties in the 



determination of the top mass (c.f. Sect.^) can 
also lead to effects of similar magnitude. 

2.2. Energy distributions 



Figure 2. Distribution of W energy for mt =174 
GeV without (dashed line) and with (solid line) 
QCD corrections 



2.2.1. Energy of W 

The dominance of the two-body t bW decay 
implies that the energy of W is strongly peaked. 
In the top quark rest frame the two-body kine- 
matics implies 



— 
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The effects of non-zero W width and QCD cor- 
rections have been studied in In Fig.|^ the 
energy distribution dr/d£'„ is shown as the func- 
tion of £'w for mt— 174 GeV. The Dirac delta spike 
in £'w is smeared by the effects of non-zero W 
width, see dashed line. QCD corrections result in 
further distortion of the spectrum (solid line) be- 
cause hard gluon radiation leads to a long radia- 
tive tail and the virtual corrections significantly 
reduce the height of the peak. 

2.2.2. Energy distributions of leptons 

The energy spectra of the leptons from semilep- 
tonic decays of the top quark are sensitive to V-A 



form of tbW vertex. In first order QCD approxi- 
mation neglecting the masses of b quark and the 
leptons one derives the following normalized en- 
ergy distributions of the charged leptons and the 
neutrinos: 



Ai (xi) 



dN 

dxf 
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dxj. J^o(y) - as-^i(y) 

(8) 

for y < Xf jy < 1, where Xi^i, = 2Ef^i,/mt de- 
note the scaled energies in the t rest frame, the 
functions J-'o{y) and J^i{y) have been defined in 
eqs.(| and (|), 



F+(x,2;) = x{l-x) (9) 
Fo-(x,y) = {x-y){l-x + y) (10) 

and [|l],|| 
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Assuming the Standard Model V-A structure of 
the charged current the spectrum of the charged 
lepton vanishes at a;^ = 1 and the spectrum of 
the neutrino does not. The latter is also signif- 
icantly harder, see solid lines in Fig.Ha-b. The 
average (scaled) energy of the neutrino (x,y) is 
larger than (xi) of the charged lepton. For the 
sake of simplicity we neglect QCD corrections and 



Figure 3. Energy distributions a) Ai(a;^) of the charged lepton and b) A^(x^) of the neutrino for the 
standard model V-A couphng (k^ ~ 0) and an admixture of V+A current (k^ =0.1) for y =0.25 and 
as =0.11. 



obtain the results shown in Table |^. It is inter- 
esting that for realistic range of nit the difference 
(x^ — xi) is significantly greater than zero. This 
observation may be useful because a value close to 
zero can be expected for the background b + W 
events of mass close to mt which do not come 
from t decays and are not correlated by the dy- 
namics. We give also the V-A prediction for the 
average difference between = 2£'> /mt of the 
more energetic and x< = 2E^/mt of the less en- 
ergetic lepton. Let us remark that for hadronic 
decays, e.g. t — > bW bdu, up-type quarks play 
the role of neutrinos and down-type antiquarks 
the role of charged leptons. Thus, (x> — a;<) 
is also proportional to the average difference be- 
tween the energies of the more and less energetic 
jets for hadronic decays of top quark. The val- 
ues given in Table || are much larger than the 
result of a crude estimation assuming isotropic 
distribution of leptons in W rest frame and sta- 
tistically independent W and b which implies 
{xy ~ x<)hckg = — v) ■ Howcver, the former 
assumption is not realistic in particular for pp col- 
liders. 



For V-f A coupling the charged lepton and the 
neutrino energy spectra would be interchanged 
in comparison to the V-A case. In |23| effects 



have been studied of a small admixture of non- 
standard V-l-A current on distributions of lep- 
tons. The thW vertex has been parametrized as 



(14) 



gv = {l + n)/^/T+^ 
gA = (-1 + k)/VTT^ 



Hence k = corresponds to pure V-A and k = oo 
to V+A. In Fig. ||a-b the lepton spectra are plot- 
ted corresponding to = 0.1, see dashed lines. 
It can be seen that the deviations from the results 
of the Standard Model (solid lines) are rather 
small. They are larger for the polarization de- 
pendent distributions which are discussed in the 
following. 

2.3. Chirality of b jet 

A measurement of b quark chirality in top de- 
cay offers another opportunity to test the V-A 
form of tbW vertex. Chirality and helicity are 
nearly identical for the highly relativistic b quark 



Table 2 



Moments of the energy distributions of the leptons in t ^ bW be^v and the Hght quark jets in 
t bW bdu decays. 
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mt=150 


mt=175 


mi=200 




1 1 All 1 ii'^ 

2(l+2y) 


0.707 


0.663 


0.631 




l+2y+3y^ 
2(l+2y) 


0.577 


0.546 


0.529 


{Xu — Xi) 


l+2y 


0.130 


0.117 


0.102 


{x> - a;<) 


3(l+2y-3y=) 
8(l+2y) 


0.317 


0.340 


0.353 



which originates from the decay, so, a high degree 
of primordial polarization is expected. 
It has been proposed long ago that distribu- 
tions of charged leptons from semileptonic decays 
of beautiful hadrons can be used in polarization 
studies for b quarks. Some advantages of neutrino 
distributions have been also pointed out p5| , p6[ . 
Recently there has been considerable progress in 
the theory of the inclusive semileptonic decays of 
heavy flavor hadrons. It has been shown that in 
the leading order of an expansion in inverse pow- 
ers of heavy quark mass I/ttlq the spectra for 
hadrons coincide with those for the decays of free 
heavy quarks ||2^ and there are no AqcD/mQ 
corrections to this result away from the energy 
endpoint. ^nn/mQ corrections have been cal- 
culated in [^,^ for B mesons and in for po- 
larized A(, baryons. For some decays the results 
are similar to those of the well-known ACCMAI 
model 1^^. Perturbative first order QCD correc- 
tions to semileptonic decays of polarized b quarks 
are also known |UJ2^,^ . An old conflict between 
|3ll and [|T| can be considered as solved in favor 
of the results of the latter work. 
The theory is therefore in good shape. The real 
drawback is that due to hadronization the net lon- 
gitudinal polarization of the decaying b quark is 
drastically decreased. In particular those 6 quarks 
become depolarized which are bound in B mesons 
both produced directly and from B* — > Bj tran- 
sitions. Only those 5's (a few percent) which frag- 
ment directly into A;, baryons retain information 



on the original polarization ]3^ , |33| | . The signal is 
significantly reduced and feasibility of polariza- 
tion studies for b jets is problematic. The on- 
going experimental studies of b polarization at 
LEP should shed some light on this problem, see 
p^ , |35t for a review. The reason is that bb 
decays can be viewed as a source of highly polar- 
ized b quarks and the polarizations of b jets in 
and t decays are similar. According to the Stan- 
dard Model the polarization of b quarks from Z'^ 
decays depends weakly on the production angle. 
The degree of longitudinal polarization is fairly 
large, amounting to {Pb) = —0.94. QCD correc- 
tions to this Born result are about 3% [^ . 

2.4. Polarized Top Quarks 

The analysis of polarized top quarks and their 
decays has recently attracted considerable atten- 
tion. Studies at a linear electron-positron collider 
are particularly clean for precision tests. How- 
ever, also hadronic |37-39| or 77 collisions |40j 
and subsequent spin analysis of top quarks might 
reveal new information. These studies will re- 
sult in determination of the top quark coupling 
to the W and Z bosons either confirming the 
predictions of the Standard Model or providing 
clues for physics beyond. The latter possibility is 
particularly intriguing for the top quark because 
mt plays an exceptional role in the fermion mass 
spectrum. 

A number of mechanisms have been suggested 
that will lead to polarized top quarks. In Bol 



this possibility has been discussed for 77 colhsions 
with circular polarized photons. Related studies 
may be performed in hadronic collisions which in 
this case, however, are based on the correlation 
between t and t decay products |^,^,^. The 
most efficient and flexible reactions producing po- 
larized top quarks are electron-positron collisions. 
A small component of polarization transverse to 
the production plane is induced by final state in- 
teractions which have been calculated in pertur- 
bative QCD | ]37| , ^ . The longitudinal polariza- 
tion Pi is large. Its dependence on the produc- 
tion angle, beam energy and the top mass has 
been discussed in . Pl varies strongly with 

the production angle, e.g. between nearly 0.6 for 
cos-d = —1 and —0.3 for cos'd = 1 at = 500 
GeV. Averaging over the production angle leads 
therefore to a significant reduction of Pl with typ- 
ical values of (P^) around -0.2 |Q. QCD correc- 
tions change (Pl) by a relative amount of about 

3% de). 

All these reactions lead to sizable polarization 
and can be used to obtain information on the 
production mechanism. However, two drawbacks 
are evident: production and decay are mixed in 
an intricate manner, and furthermore the degree 
of polarization is relatively small and depends on 
the production angle. Top quark production with 
longitudinally polarized electron beams and close 
to threshold provides one important exception: 
the restricted phase space leads to an amplitude 
which is dominantly S-wave such that the elec- 
tron (and positron) spin is directly transferred 
to the top quark. Close to threshold and with 
longitudinally polarized electrons one may deal 
with a highly polarized sample of top quarks in- 
dependent of the production dynamics. Thus one 
can study t decays under particularly convenient 
conditions: large event rates, well identified rest- 
frame of the top quark, and large degree of polar- 
ization. 

2.4.1. Angular distributions 

In the rest frame of the decaying t quark the an- 
gular distributions of the decay products are sen- 
sitive to its polarization. Let us define the angle 
Ov, between W boson three-momentum and the 
polarization three-vector s. In the top quark rest 



frame s = (0, s ) is the polarization four-vector of 
the decaying top quark. Note that S — \s \ — 1 
corresponds to fully polarized and S' = to un- 
polarized top quarks. We define also the angles 
Qj^ and between s and the directions of the 
charged lepton and the neutrino, respectively, and 
0< for the less energetic lepton in semileptonic or 
less energetic light quark in hadronic decays. For 
the sake of simplicity let us confine our discussion 
to Born approximation and consider semileptonic 
t bW M^iy and hadronic t bW hdu 
decays. 

The angular distribution of the charged lepton is 
of the form 



a cos 0+ 2 



(15) 



which follows, c.f. the following subsection, from 
factorization of the angular-energy distribution 
into an energy and an angular dependent part. 
This factorization holds for arbitrary top mass 
below and above the threshold for decays into 
real W bosons [43 2^. It is noteworthy that for 
S=l the angular dependence in (^) is maximal 
because a larger coeffecient multiplying cos 0+ 
would be in conflict with positivity of the decay 
rate. Thus the polarization analysing power of 
the charged lepton angular distribution is maxi- 
mal and hence far superior to other distributions 
discussed in the following. In particular the an- 
gular distribution of the neutrino reads |2^ ]: 

^4[l + M.)5cos.,] (16) 

where hy(y) is given in Table |^. The distribu- 
tion of the direction of W can be easily obtained. 
Only the amplitudes for the helicity states of W 
Aw = ~1 and A„ = are allowed and their contri- 
butions to the decay rate are in the ratio 2y : 1 
Q] . The corresponding angular distributions are 
of the form 



diV-1,0 ^ 1 
d cos 6„ 2 



(1t5cos6'w) 



(17) 



After summation over the W polarizations the 
following angular dependence is obtained: 



dN 



1 



dcos( 



[1 -I- h„{y)S cos 9y, 
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Table 3 



Angular depenaence or tne distributions or W bosons, neutrinos and less 
be~^iy or light quark jets in i ^ bW — > bdu decays. 


energetic leptons in t 


bW 


y 


TOt = 150 


mt—175 


mi=200 


U 1 12y{l~y+y\ny) 

I^Ay) ^ (l-a)^(i+2y) 


-0.521 


-0.311 


-0.127 


"-wyyj l+2y 


0.275 


0.410 


0.515 


U 1 6y{l-y-2yln[{l+y)/{2y)]} 


0.464 




n 'I'lQ 



where h^{y) is also given in Table I It IS evi- 
dent that the charged lepton angular distribution 
is significantly more sensitive towards the polar- 
ization of t than the angular distributions of W 
and V. The charged lepton is likely to be the less 
energetic lepton because its energy spectrum is 
softer than that of the neutrino. For large values 
of rrit the angular distribution of the less energetic 
lepton 



d cos 2 



(19) 



is a more efficient analyser of top polarization 
than the angular distribution of neutrinos. For 
rrit in the range 150-200 GeV it is also better than 
the direction of c.f. Table |. 

2.4.2. Angular-energy distributions 

For semileptonic decays the normalized distri- 
butions of leptons including first order QCD cor- 
rections can be cast into the following form: 

, = \ [Ai(x,) + 5cos0+Bi(a;^)] 

ax I d cos t'-|_ 2 

(20) 

- — = \ [K{xu) + S'cos6'oB^(a;^)] 
axv a cos Oq 2 

(21) 

where the functions A\{x) and A^(x) have been 
defined in eqs. (0) and (|), and 

„ , > 'fn(x,v) — 37 (x,y) 

B.x = 12 ^^^^^^ -rrv^ 23 



where 

^o{x,y) -- 
and [Esira 



^t{x,y) (24) 
{x-y){l-x + y-2ylx) (25) 



3t{x, y) = J([(a;, y) $o - x<^i + - 2x - y 
-2y/x - 2/2;)$2o3 + (-3 + x -f 2/x)^i 
+ {-2xy + 3x~ Ax'^ + 6y - y'^ - 2y^/a;)$5 
+ (2 - 2x^ + 2y- 3y^ - 2y/x + 3y^/x)/2 

(26) 

Ji{x,y) J^(x,2/) ^0 + {-2xy - X - 5y+ 
- 2y^/x)'i>i + (3 + 10a; + y + lOy/x 
-2/a;)$2o3 + (-3 + 12a;y + x - 7y - y^ 
-12y/x + 8y'^/x + 2/x)$4 + (62;?; - 9x 
-4x^ -y-2y^ + 10y^/x)^5 + (2 - 5xy 
-2x^ + 2y + 7y^ - 2y/x-2y^/x)/2 (27) 

Eq.(p^ implies that in Born approximation the 
double differential angular-energy distribution of 
the charged lepton is the product of the energy 
distribution and the angular distribution (|l5|). 
QCD corrections essentially do not spoil this fac- 
torization For the neutrino such factoriza- 
tion does not hold, c.f. eqs.([lo|) and (|25|). After 
integration over the angular dependence of the 
neutrino distribution is much weaker than for the 
charged lepton, c.f. eq.(|l^). 

In Fig.||the functions Bi(x) and Bi,{x) are shown 
as solid lines for y —0.25 and =0.11 |^^. The 
effect of non-standard coupling defined in eq.([l^) 
is much stronger for the polarization dependent 



Figure 4. Angular-energy distribution functions in the Standard Model (k^ — 0) and for the admixture 
of V+A current (k^ =0.1): a) ^\(xi) for the charged lepton and b) ^i,(xy) for the neutrino, j/=0.25 and 
a^=0.11. 



distribution of neutrinos, see dashed lines in Fig.Q 
corresponding to = 0.1 In Table ^ the 

moments 

Jy 

Sf") = / d2/x^-B,,,(x,y) (28) 
Jy 

are given for integer k between -1 and 3, y = 0.25 
and as{mt) = 0.11 Q . The upper entries in 
the table denote the values of the moments for 
= and the lower ones the ratios of the mo- 
ments evaluated for = 0.1 to those for = 0. 
It is evident that the moments b]^^ which gov- 
ern the angular dependence of the neutrino spec- 
trum are particularly sensitive towards a V-f A 
admixture. The effect is most pronounced for 
the moment fc = — 1 which enhaces the lower en- 
ergy part of the spectrum and where the relative 
change amounts to 17% for = 0.1 . Thus, the 
angular-energy distribution of neutrinos from the 
polarised top quark decay will allow for a partic- 
ularly sensitive test of the V-A structure of the 
charged current. The effect of QCD correction 



can mimic a small admixture of V+A interaction. 
Therefore, inclusion of the radiative QCD correc- 
tion to the decay distributions is necessary for a 
quantitative study. 

3. TOP QUARK PAIR PRODUCTION 

3.1. e+e tt 

It is evident that the bulk of top studies at 
an e^e~ collider will rely on quarks produced in 
e+e~ annihilation through the virtual 7 and Z, 
with a production cross section of the order of 
CTpomtQ- For quarks tagged at an angle i?, the 
differential cross section in Born approximation 
is a binomial in cos -d 

= 1(1 + cos^ d)au + -, sin^ ^ul 
dcosv 8 ' 4 

3 

-I- - cos zJcTi. (29) 

\J and L denote the contributions of unpolarized 
and longitudinally polarized gauge bosons along 
the d axis, and F denotes the difference between 

^This subsectionis a shortened version of a comprehensive 
review given in |12|. 



Table 4 

The moments A^^''^"^ and B^^''^'^ of the angular-energy distributions for the charged leptons and for the 
neutrinos for y = 0.25 and as = 0.11 for — (upper entries) and the ratios between the moments for 



= 0.1 and 


= flower entries). 








k 










-1 


2.008 
.981 


2.005 
.940 


1.593 
1.023 


-0.707 
1.166 





1.000 
1.000 


.998 
.949 


1.000 
1.000 


-0.452 
1.110 


1 


.559 
1.021 


.558 
.960 


.683 
.984 


-0.322 
1.068 


2 


.345 

1 0/1'^ 
1.U40 


.344 


.500 
.y / o 


-0.249 

1 0*^7 
1 .Uo ( 


3 


.230 
1.064 


.230 
.989 


.385 
.965 


-0.203 
1.015 



right and left polarizations. The total cross sec- 
tion is the sum of U and L: 



(30) 



In Born approximation the coefRcients ct* can be 
expressed in terms of the cross sections for the 
massless case 



(Jq = per + p a 



'B 

where 



/3= v/l-4m2/s 
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and 

D (s, m^, r^) = (s - + (sTz/toz)^ 



The fermion couplings are given by 



2ll - 4e/ sin^ 9. 



a/ 



21: 



f 



(35) 
(36) 
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and the possibihty of longitudinal electron po- 
larization {p — — 1;-|-1;0 for righthanded; left- 
handed; unpolarized electrons) has been included. 
QCD corrections to this formula are available for 
arbitrary m^/s up to first order in as- 



(3~/3^).,y 
2 



Ev-K/JV-^^ Ra 



(38) 



The exact result [Q for the QCD enhancement 
factors can be well approximated by 



Rv 



1 



3 



TT 

2p 



(3 ^TT 



3 

in 



(39) 



Figure 5. Cross section for tt production in units of Upoint including resonance and QCD enhancement 
and initial state radiation. 



and 

Ra = 

Aas r TT /19 22 ^ T r.2\ f-^ 3 \" 

i + -b-(io-T^+2^) (2-4^) 

(40) 

For small (3 these factors develop the familiar 
Couloumb enhancement ~ '^iras/fi compensat- 
ing the phase space suppression ~ /3. This leads 
to a nonvanishing cross section which smoothly 
joins the threshold region. 

Initial state radiation has an important influence 
on the magnitude of the cross section which is 
significantly suppressed in the threshold region. 
The correction factor increases rapidly with en- 
ergy, but stays below 0.9 in the full range un- 
der consideration. In Fig.|^ the cross section for 
e+e^ — > ti is shown as the function of the center- 
of-mass energy Q. 

3.2. Pair production near energy threshold 

It follows from Table |l| that the top quark is a 
short-lived particle. Its width Tt is in the range 
1-2 GeV and fairly exceeds the tiny hyperfine 
splitting for open top hadrons, the hadronization 
scale of about 200 MeV, and even the energy split- 
ting between 15* and 2S tt resonances. On one 
side this is an advantage because long-distance 
phenomena related to confinement are apparently 
less important for top quarks pst . In particular 



depolarization due to hadronization is practically 
absent. On the other side the amount of informa- 
tion available from the threshold region is signif- 
icantly reduced. Toponium resonances including 
15 overlap each other. As a consequence the cross 
section for ti pair production near energy thresh- 
old has a rather simple and smooth shape. Nev- 
ertheless, as first pointed out by Fadin and Khoze 
the excitation curve a{e^e~ tt) allows a 
precise determination of rrit as well as of other 
physical quantities such as Tt and the strong 
coupling constant. The results of for the 
Coulomb chromostatic potential have been de- 
rived analytically in the non-relativistic approxi- 
mation. Strassler and Peskin have obtained 
similar results using numerical approach and a 
more realistic QCD potential. The idea p9| , ^ to 
use Green function instead of summing over over- 
lapping resonances has been also applied in cal- 
culations of differential cross sections [^-^ . In- 
dependent numerical approaches have been dev- 
oleped for solving Schrodinger equation in the po- 
sition space [0 and Lippmann-Schwinger equa- 
tion in the momentum space |5^,|5^. The results 
of these two methods agree very well [Q. One 
of the most important future applications will be 
in determination of mt and as- It has been ar- 
gued js^] and demonstrated [^,^ that the com- 
bined measurements of the total and the differ- 
ential cross sections in e^e~ — > tt offer a very 



Figure 6. Distribution \pG{p, E)\'^oi the momentum p — \p\ oi the top quark in tt system for E=-2.9 GeV 
(15' peak), E=0 and E^2 GeV a) mt= 150 and b) 180 GeV. 



promising method for a simultaneous determina- 
tion of TOt and as{mt). The estimated errors of 
such a determination are js^: Amt—300 MeV 
and Aas=0.006 for mt==150 GeV and ATOt=520 
MeV and Aas=0.009 for mt=180 GeV. The er- 
rors are correlated, so for as fixed by other mea- 
surements a much better determination of 
can be obtained. In the following the method of 
|52|j5^ ] is presented. Some theoretical uncertain- 
ties related to the momentum dependent width of 
tt system, long-distance part of the QCD poten- 
tial and the definition of the pole top quark mass 
are discussed. Efi'ects of Higgs exchange are also 
considered. 

3.2.1. Green function 

Let us describe briefly the Green function 
method for e'^e~ — > tt annihilation near the 
energy threshold and the numerical solution of 
the Lippmann-Schwinger equation; see p2| , p3[ 
for details. We simplify the original calculations 
and neglect the Z'^ contribution and the trans- 
verse gluon correction to the production vertex. 
The differential cross section for the top quark 
pair production in electron-positron annihilation 
reads: 



da 
(Pp 



(p, E) - ^^r(p, E) \g{p, Ef (41) 



where r(p, E) — Ttt/2 . The Green function 
g {p, E) is the solution of the non-relativistic 
Lippmann-Schwinger equation 



g{p,E)^go{p,E) 

Qoip.E) 
where 



^,viP~?)gim 



E = ^-2mt 



(42) 



(43) 



denotes the total energy of the ti system and p 
the momentum of t quark. V{p) is the potential 
in momentum space and T{p, E) — T^i/2, where 
denotes the width of the the ti system. The 
free Hamiltonian that is used to define the Green 
function g^ includes the momentum dependent 
width: 

goip, E) = (44) 

E~^^+iTip,E) 

Near the energy threshold one can neglect all but 
S partial waves and solve numerically the corre- 
sponding one-dimensional integral equation. The 
spherically symmetric solution fulfills the unitar- 
ity condition | ^ 

d^p 



Tip,E) \gip,E)\' = 
-ImG{x = 0,x' = 0,E) 



(45) 



which for the constant decay rate reduces to the 
formula for the total cross section derived in 
P^ , ^ . Eq.(^ can be also considered as a non- 
trivial cross check of the numerical approach. 
In Fig.^ the momentum distributions of t quarks 
in ti systems of energy E corresponding to 15 
peak , E=0 and E=2 GeV are shown as solid, 
dotted and dashed lines respectively. The shift 
towards larger momenta with increasing E and 
the narrowing of the distributions is clearly visi- 
ble. It has been shown ||5^,|5^ that the position of 
the maximum is not sensitive to the initial state 
radiation, so it can be used to measure E and 
rrit- On the other hand the rapid increase of the 
cross section a{e~^e~ ti) is related to the lo- 
cation of IS' resonance which depends on both 
rut and as- Thus, the combined measurement of 
both the total and differential cross sections leads 
to a simultaneous determination of these param- 
eters. The results shown in Fig.^ have been ob- 
tained in [^2| assuming constant width of ti sys- 
tem r(p, E) —Vf Let us discuss now theoretical 
problems which appear when this assumption is 
lifted. 

3.2.2. Width of ti system. 

The width of the ti system depends on the 
momentum of t quark because both the matrix 
element and the phase space available for the 
decay products depend on it. When produced 
near energy threshold t and i cannot be consid- 
ered as free particles. The binding energy and 
the kinetic energy of internal motion tend to re- 
duce the available phase space for the decay. Al- 
though the effect is only 0{as^) the suppression 
can be significant | |56[ |. Thus, in a high preci- 
sion calculation one has to consider the width 
T^_i{j>) as a non-trivial function of the momen- 
tum p. The phase space effect tends to reduce 
the decay rate of bound top quarks relative to free 
ones [ ^6|j5l|] and the effect is enhanced because for 
short-lived particles the momentum distribution 
is broad. However, for the same reason the de- 
cays take place at short relative distances, where 
the wave functions of b and h quarks originat- 
ing from the decays are distorted {enhanced) by 
Coulomb attraction. Therefore, when calculating 
the amplitude of i — > hW transition, one should 



Figure 7. Comparison of the annihilation cross 
section (^(e+e^ ti) evaluated for the constant 
width Tt (dashed line) and the momentum depen- 
dent widths including time dilatation (solid line) 
and phase space suppression (dotted line). 



use Coulomb wave functions rather than plane 
waves for h quarks. This effect clearly increases 
the rate. A third factor is due to time dilatation: 
a top quark moving with velocity v lives longer in 
the center-of-mass frame. While phase space re- 
duction and time dilatation can be implemented 
in a straightforward way Coulomb enhancement 
cannot be easily taken into account. In principle 
one has to replace the plane wave functions for 
b quarks by relativistic Coulomb functions when 
evaluating the amplitude for the t bW tran- 
sition. In Figj^ |^3| predictions for the annihila- 
tion cross section a^e'^e^ —^ ti) are shown corre- 
sponding to different asumptions about Ttf. solid 
line has been obtained including only time dilata- 
tion, dotted line including only the phase space 
suppression and the dashed line corresponds to 
the constant width. The effect of the phase space 
suppression is non-negligible whereas time dilata- 
tion produces rather small effect. 
It has been conjectured |]5^ , p3| , ^ that, in close 
analogy to what happens in the case of negative 
muons bound in nuclei for chromostatic at- 
traction in ti systems the phase space suppression 
and the Coulomb enhancement cancel each other 
and only the suppression due to time dilatation 
should be included. Recently this conjecture has 



been proven |59|| . Moreover, it has been explicitly 
demonstrated that unphysical gauge dependence 
of the momentum dependent width [Sl]] also dis- 
appeares in the final result ||5^. The distribu- 
tions of the top quark momentum, see Fig.^, are 
also only weakly affected. It should be stressed, 
however, that the experimentally accessible three- 
momentum distributions of Wb systems from top 
quark decays are not identical to the distributions 
of the top quark momentum. Calculable correc- 
tions arise due to final state rescattering ^ . 

3.2.3. QCD potential 



Figure 9. QCD potential VjKxi''') for different 
values of asim^): solid - 0.12, dashed - 0.11, 
dash-dotted - 0.13 and dotted - 0.10 and 0.14. 



Figure 8. aeff{q) for different values of Q;s(mz): 
solid - 0.12, dashed - 0.11, dash-dotted - 0.13 
and dotted - 0.10 and 0.14. 



Some difficulties appear also for the QCD po- 
tential. In order to illustrate these problems let 
us consider the following potential defined in 
momentum space {p = |p |): 

IQlTOeffip) 



Vjkt{p) 
The function 
aeffip) = 



5 

apert (P) 



Vo 6{p) 



if P > Pi 
if p < P2 



(46) 



(47) 



is shown in Fig.^. At large momenta, i.e. for 
p > pi—5 GeV Vjkt{p) is equal to the pertur- 
bative two-loop QCD potential |60) for Uf = b 
quark flavours and as{mz)= 0.10-0.14. At inter- 
mediate and small momenta, i.e. for p < p2=2 
CeV a Richardson-like phenomenological poten- 
tial ||6l| is employed, and a linear interpolation 
formula is used in between. While the form of 
the potential in momentum space is fixed at large 
and intermediate momenta by perturbative QCD 
and phenomenology of bb and cc quarkonia, it is 
not known in the infrared (confinement) region. 
Thus Vq in eq.(46) is to some extent an arbitrary 
parameter. In |56| Vq has been fixed by imposing 
an evidently arbitrary condition on the potential 
in the position space 



VjKT{r = IGeV"^) = -l/4GeV 



(48) 



Then after Fourier transformation the potential 
VjKT{r) is obtained, see Fig.^, which describes 
well the bb and cc quarkonia. A different choice of 
Vq results in a redefinition of b and c quark masses 
which in the framework of the non-relativistic ap- 
proach are defined only with limited (~ 300 MeV) 
precision. In particular one can impose conditions 
relating Vq and as{mz) which lead to different en- 
ergy dependence of the momentum distributions 



in systems. At this point a brief discussion is in 
order on the relation between the top masses de- 
termined on the basis of different potentials j62) . 
The perturbative two-loop QCD potential in mo- 
mentum space is fixed unambiguously for suffi- 
ciently large Q^. In order to calculate the po- 
tential in coordinate space the small behavior 
has to be specified in an ad hoc manner. Dif- 
ferent assumptions will lead to the same short 
distance behavior. The potentials will, however, 
differ with respect to their long distance behav- 
ior. In it has been argued convincingly that 
the long distance tail is cut off by the large top 
width. However, an additive constant in coordi- 
nate space can be induced by the small momen- 
tum part of V{p). This additional term leads to 
a shift in the tt threshold, which in turn can be 
reabsorbed by a corresponding shift in rrit- The 
different assumptions on the long distance behav- 
ior are reflected in differences between the predic- 
tions of 1 50 5^j5^ for the precise location of the tt 
threshold for identical values of and as well 
as in differences in the dependence of the mo- 
mentum distributions for fixed rrit and energy, see 
also |54|. All these differences can be attributed 
to the freedom in the additive constant discussed 
before. The same additive constant appears in 
bb spectroscopy, such that the mass difference be- 
tween top and bottom is independent from these 
considerations. Recently problems with the defin- 
ion of b quark mass due to long distance effects 
(infrared renormalons) have been noted in[Q. 

A similar study has been performed in ||64|. In 
Fig.|l^ (T(e+e~ — > ti) is shown as the function 
of energy relative to the location of the 15 peak 
AE' = ^/s — Els ■ The solid lines have been ob- 
tained for the potential of ||5^ and the dashed 
ones for the potential of ||5^. These two poten- 
tials are equal at short distances but different pre- 
scriptions are used at intermediate and large dis- 
tances. It can be seen that the height of the IS* 
peak is affected. The ambiguity may be fixed 
by performing a simultaneous fit to the differ- 
ential momentum distribution, angular distribu- 
tions (e.g. the forward-backward asymmetries) 
and/or by fixing as from independent measure- 
ments. It is clear that this problem deserves a 
careful further study. 



Figure 10. Comparison of tt cross section as the 
function of energy relative to the position of 15 
peak for different potentials |64|. 



3.2.4. Higgs effects in the threshold region 

Effects of Higgs exchange may be non- 
negligible in the threshold region |]50| in particular 
if the Higgs mass is relatively low as predicted by 
supersymetric models. The dominant effect from 
a light Higgs boson can be described by an in- 
stantaneous Yukawa potential 



Vvukir) = -K- 



with 



\/2Gi=^TO?/47r 



(49) 



(50) 



The vertex correction from Higgs exchange can 
be included in the potential calculation by adding 
Vyuk to the QCD potential: 



^e// = VqcD + Vyi 



uk 



(51) 



Then, Veff can be used instead of Vqcd in nu- 
merical calculations of the total and differ- 
ential ||5l|,|5^ cross sections for e+e~ annihilation 
near tt threshold. The dominant contribution to 
the remainder can be included as an overall cor- 
rection factor to the amplitude pq]: 



1 



fthr{ml/ml) 



ruh 



(52) 



where |66h6 



fthrir) = -— [-12 + 4r+(-12 + 9r-2r2)lnr 



+ -(-6 + 5r-2r^)Ur) 
r 



r arccos 



(V^/2) 



(53) 



if r < 4 



-Jr{r~A)\\u '+^lEIE ifr>4 

(54) 



This factor is energy independent and can be used 
below as well as above the threshold. 
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